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Here, we report results of an experiment creating a transient, highly correlated carbon state using a 
combination of optical and x-ray lasers. Scattered x-rays reveal a highly ordered state with an electrostatic 
energy significantly exceeding the thermal energy of the ions. Strong Coulomb forces are predicted to induce 
nucleation into a crystalline ion structure within a few picoseconds. However, we observe no evidence of 
such phase transition after several tens of picoseconds but strong indications for an over-correlated fluid 
state. The experiment suggests a much slower nucleation and points to an intermediate glassy state where 
the ions are frozen close to their original positions in the fluid. 



Carbon is known to have many phases with extraordinary properties including the highest melting tem- 
perature of any material. Yet, the high-pressure transitions into phases with temperature above 5,000 K 
and densities above 1 g/cm^ are poorly understood^'^. Carbon is expected to form a chemically inert phase 
under such extreme conditions, and exists in liquid form in the interiors of giant planets, such as Neptune, 
Uranus, and carbon-rich extrasolar planets^'^. The high-pressure phases of carbon also determine the state and 
evolution of more compact objects like white dwarf and neutron stars^'^. 

Strongly driven systems far from equilibrium offer a very direct way of investigating the microscopic dynamics, 
inter-particle interactions and phase changes. High-pressure melting^, transient changes of optical properties''^ 
and the relaxation path towards equilibrium^^'^^ are examples with high importance when studying materials with 
high energy density. Indeed, such states that naturally occur in large astrophysical bodies can only be created in the 
laboratory by depositing a large amount of energy into a sample and driving it far from equilibrium. However, the 
study of such warm dense matter (WDM) states poses severe challenges in both modeling and experimental 
diagnostics. Their temperature is too high for solid state theories to be applicable, but at the same time WDM is too 
cold for plasma perturbation techniques to converge since the thermal energy is only a fraction of the Coulomb 
energy between the particles^^. Accessing WDM states experimentally has also remained rather challenging due to 
the inability to directly probe stellar matter and, in the laboratory, due to the difficulty in achieving well-diagnosed 
off-Hugoniot states created by large lasers and pulsed-power machines^^ While some experiments have produced 
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a transient liquid state of carbon through isochoric heating^, only 
limited structural and dynamical information could be obtained. 

To overcome these limitations, we used here two nearly collinear 
laser beams to create a unique dense matter regime in an experiment 
performed at the Linac Coherent Light Source^^. First we illuminated 
a 1 jim thin graphite foil with an optical laser (see Figure 1 for details 
of the experimental setup). The laser was focused with a relatively 
large focal spot of 80 jim diameter and a moderate intensity of 
— 10^5 W/cm^. The corresponding ablation pressure of —30 Mbar^^ 
launches a shock wave through the carbon sample. The shock com- 
presses the sample to densities of p ~ 2.5-5 g/cm^ and heats it to 
temperatures T ~ 5, 000-10, 000 K, inducing the carbon to melt (see 
also Ref. [2] for the high pressure phase diagram of carbon). Figure 2 
shows an example for the density and temperature profiles after 
40 ps as obtained from modelling the radiation driven hydrodyn- 
amics. After a time of ^ ~ 40-50 ps, the shock reaches the opposite 
side of the foil and breaks out. This long wavelength, optical laser also 
produces a large number of nonthermal electrons having a high- 
energy tail with a temperature of —7.5-14 keV^^. These electrons 
will significantly enhance the ionization degree of the carbon sample 
by electron impact. 



After the illumination by the optical laser, the x-ray beam from the 
free electron laser (FEL) interacts with the pre-formed liquid state at 
different delays. The high intensity and photon energy of the FEL, 
— 10^^ W/cm^ and E^ = S keV, leads to additional ionisation of the 
sample^'^. While the FEL intensity is only sufficient to photo-ionize 
about 0.1% of the ions in the probed volume directly, the energy of 
the ejected electron, being the difference between the photon energy 
and the K-shell binding energy of the ions, is large allowing for the 
creation of several secondary electrons contributing roughly 10% of 
the total ionization of the sample. 

Photo and impact ionisation are enhanced by Auger recombina- 
tion of core holes that are filled within a few femto- seconds and 
results in the ejection of an L-shell electron. Therefore, the ionization 
of the sample is expected to be far above the charge state of a system 
in thermodynamic equilibrium. Given that electrons with energy 
>7.5keV will perform many ionization events during the FEL 
pulse^^ and have a range of less than 1 jim, most of the carbon sample 
in the probed region will consist of highly charged ions. The hot 
electrons also deposit additional energy into the background plasma, 
raising the electron temperature by a few eV in a time shorter than 
the duration of the FEL pulse. On the other hand, electron-ion 
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Figure 1 | Experimental setup. This work was been performed at at the LCLS XPP instrument. A 300 fs long, 10 mj TiiSapphire laser (operating at 
wavelength of 800 nm) was focused onto a 1 [im thick graphite foil, mounted on a rastering stage within a vacuum chamber, to a focal spot-size of 80 |am 
diameter resulting in an intensity on target of 10^^ W/cm^. The shocked heated carbon foil was then illuminated by the FEL probe beam (operating at 
8 keV, with ~ 1 mJ energy in a 80 fs long pulse) propagating nearly coUinearly with the optical laser pulse. The x-ray spot size was around 20 [im 
diameter, focused using a series of beryllium refractive lenses. The inset shows the spectrum of the x-ray beam as it reaches the sample, as well as the 
measured scattering spectrum atO = 50°. Temporal synchronization and spatial overlap of the optical and x ray beams was achieved through the use of 
optical damage shadowgraphy^^. The scattered x-rays are collected using a highly oriented pyrolytic graphite (HOPG) crystal spectrometer in von Hamos 
geometry. Details of the crystal spectrometer calibration are given in Ref. [23] . Data is then recorded onto a high quantum efficiency, high repetition rate 
pixel-array detector (CSPAD). An example of the raw data recorded on the detector is shown in the figure. The diffracted photon energy increases from 
bottom to top in the image. From the raw image, the spectrum was obtained by integrating along the non-dispersive direction, but only in a narrow central 
strip, as indicated in the figure. This is to minimize spectral broadening associated to crystal aberrations, as discussed in Ref. [23]. The spectrometer, 
including detector, was mounted to a six- axis robotic arm, enabling the angle at which the scattered radiation was measured to be varied during the 
experiment (6 = 20°, 35°, 50°, 130°). The HOPG crystal was large enough to collect rays from ±5° of the nominal scattering angle. The polarization of the 
x-ray beam was at 90° to the angle of scatter, maximizing the scattering efficiency and obviating any polarization effects on the scattered radiation. 
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Figure 2 | Radiation hydrodynamic simulations. Numerical simulations of the optical laser interaction with the solid carbon foil are performed using the 
code NYM^^. The simulations have been done in cylindrical 2D symmetry, using an inverse bremsstrahlung model for laser absorption and flux-limited 
diffusion for hot-electron transport, nym simulations do not include the effect of hot electron production at the laser spot nor the FEL heating. Panel a: 
contour plot of the mass density at t = 40 ps; Panel b: electron temperature at t = 40 ps after optical laser arrival time. 



equilibration occurs on much longer times - typically hundreds of 
picoseconds - and thus the carbon ions remain at the temperature 
created by the Shockwave, that is, ~ 5, 000-10, 000 K. 

Results 

The FEL beam delivers enough photons that a sufficient number are 
scattered by the electrons in the plasma and can be used as diagnost- 
ics^^'^°. The power of the scattered photons is given by dP/Po = 
n/€(77^tot(^) where Pq is the incident power of x-rays, Ui is the 
ion density, Gt is the Thomson cross section, dQ is the acceptance 
solid angle of the analyzer crystal (see Figure 1), and Stot(^) is the total 
electron structure factor. Stot(A;) reflects the spatial arrangements of 
the electrons in the sample and, thus, relates the scattered power 
spectrum to the material properties. It is probed at a wavenumber 
oik = {47i/X) sin(^/2), where X = 1.5 A (8 keV) is the x-ray wave- 
length and 0 is the scattering angle. 

As plotted in the inset of Figure 1, the scattered radiation shows 
two distinct features: i) an elastic (or nearly elastic) peak at the 
original wavelength, and ii) a frequency- shifted feature. The first 
feature is the result of scattering by electrons that kinematically 
follow the ions^^'^^, which in WDM can be bound or form the screen- 
ing cloud around the ions. The second feature, which is smaller in 
intensity but still contributing a significant amount of the total signal 
(see Figure 1), results from inelastic scattering off free electrons and 
Compton scattering from core electrons which causes photo-ioniza- 
tion. We used a highly oriented pyrolitic graphite (HOPG) crystal 
analyzer^^ in order to separate the elastic from the inelastic scattering. 
The measured spectra were furthermore deconvolved by taking into 
account the measured spectral profile of the FEL source. Any diffuse 
scattering arising from the polycrystalline structure of the cold (un- 
shocked) region of sample is removed by subtracting the measured 
scattering signal without optical laser heating (i.e., FEL only) 
weighted by the mass fraction of unshocked material. At any given 
time the latter is estimated from the results of a radiation-hydrodyn- 
amics simulation (see Figure 2). The cold scattering signal is plotted 
in Figure 3a. The cold graphite structure is hexagonal close-packed 
(hep). The (012) diffraction line contributes to the measured peak at 
3.4 A"^ (6 = 50°). 

Information about the ion structure in the sample is contained in 
the coherent nearly elastic scattering peak. The corresponding part of 
the electron structure factor is given by \f{k) + q{k)\^Sii{k), where f{k) 



refers to the atomic form factor of bound electrons and q{k) to the 
screening cloud of electrons around an ion^\ In our experiment, the 
sum of/(A;) and q{k) is only weakly dependent on the details of the 
plasma conditions as we have klxF ~ 1 (>lrF is the Thomas -Fermi 
screening length) and kais ^ 1 (ai^ is the radius of the Is electrons in 
carbon). Thus, q{k) is very small and/(A;) ~/(0) is almost constant. 
The main modulation of the elastic scattering signal results from the 
ion-ion structure factor, Sii(A;), which describes the ion positions in 
Fourier space. To account for the detector response and transmission 
efficiency of the x-ray collection optics, we have used as the normal- 
ization condition that for A: > 7 A" ^ the structure factor should tend 
to unity. This diffraction technique, recently implemented to study 
correlations in dense aluminum^^, is different from the commonly 
used angle-resolved, frequency- integrated methods. The latter are, 
however, unable to provide a direct measurement of the ion structure 
if inelastic scattering is large, as shown in Figure 1. 

Our results for the ion-ion structure factor extracted from the 
elastic scattering feature are shown in Figure 3b, with the error ana- 
lysis presented in Figure 3c, for various delays between the optical 
laser and the x-ray probe. The measured rise of the structure factor at 
small k, as well as a finite value at large scattering angles is typical for a 
strongly coupled liquid. The presence of a high correlation peak 
indicates that the carbon ions exhibit a pronounced short-range 
order on the scale of 2n/ko ~ 2 A with ko ~ 3.4 A"^ being the position 
of the first peak in the structure factor for the conditions in our 
experiment. While the cold data also show a peak at the same posi- 
tion, the fact that the diffraction intensity is much stronger when the 
sample is heated, and the appearance of non zero scattering at —7.4 
A"^ (^ = 130°) both suggest a substantial change from the initial hep 
carbon structure. This new order is a result of the strong Coulomb 
repulsion between the highly charged ions. Indeed, such a structure 
factor is indicative for systems where the Coulomb energy between 
the ions significantly exceeds their thermal excitation. The ratio of 
these two energies is usually expressed by the parameter r = 
/ 47ieQakBT:^l with a = (3/47i/i/)^^^ being the mean distance 
between ions (where Ui is the ion density). Moreover, the total pres- 
sure of such systems, expressed as = p^^ + pex^ is dominated by the 
excess pressure due to ion-ion repulsion and not by the ideal contri- 
bution, pid, as in a gas or typical plasma^^ The Coulomb forces 
between ions thus dominate the thermodynamics of the system. 
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Figure 3 | Structure factor of strongly coupled carbon. Panel a: scattering data from the un-shocked (cold) graphite. The position and relative intensity 
of the bcc lattice peaks, labelled according their Miller indices, for p = 2.5 g/cm^ (blue lines) and p = 4.5 g/cm^ (orange lines) is indicated in the figure. 
The Bragg peak intensity is estimated by assuming rjY = 0.1. Panel b: data were collected at a repetition rate of approximately a shot per 10 seconds, 
allowing for the coherent scattering to be taken at various combinations of beam delays and measurement angle. Each data point consists of an average of 
about 50-60 individual shots. The FEL energy has been individually recorded for each shot, and then single spectra have been weighted by the 
corresponding energy in the FEL beam (with cold scattering removed). The experimental values for Saik) were derived from the total elastic scattering by 
assuming Z = 4.5 ± 0.5. The uncertainty in Z determines the vertical error bars. The standard deviation in signal intensity from shot-to-shot variations in 
a data set is within the reported errors. Horizontal error bar are related to the finite acceptance angle of the spectrometer. The calculated curves are results 
from a fluid model for charged ions embedded in a polarisable background of electrons for a density of p = 2.5 g/cm^ and p = 4.5 g/cm^, an ion charge 
state of Z =4.5, and two different temperatures ( 1 eV equals 1 1 ,600 K) . Panel c: plot of = | - Sf^^ \ ^ versus Tand Z, where S^"^ = P""^ / [f (k) + q{k)\^ 
is the structure factor derived from the measured data (P^'p) and S^^^ is the calculation for p = 2.5 g/cm^, 6 = 50° and 40 ps delay. 



In the case of a plasma in thermodynamic equilibrium, density 
functional theory coupled to molecular dynamics (DFT-MD) simu- 
lations represent an ab initio method to calculate the ion-ion struc- 
ture factor^^. The result of such simulations is shown in Figure 4. We 
notice that for the conditions of our experiment, DFT-MD gives an 
ion-ion structure factor which is relatively flat in the region of 
k> 3 A"^ indicating a moderately coupled system only. By compar- 
ing the DFT-MD results with the experimental data (Figure 3), we see 
that the agreement is rather poor, particularly in the region of the first 
peak of the structure factor. This is due to the fact that an eqmlibrium 
calculation only predicts a moderate ionization for carbon (Z ~ 2), 
whereas in our experiment, additional ionization is driven by supra- 
thermal electrons whfle keeping the initial density and ion temper- 
ature almost unaffected. Hence, the plasma state created in our 
experiment is better described in terms of positively charged ions 
with a much higher charge state embedded in a polarisable back- 
ground of electrons^^. This situation can often be described by the 
mean spherical approximation (MSA) making the problem of the ion 
structure analytically solvable. The best fit with the experimental data 
is obtained for an ion charge state of Z = 4.5 (Figure 3). The calcula- 
tions shown in Figure 4 clearly indicate that as the ionization 
increases there is the emergence of a strong correlation peak. The 
position of the peak depends on the density, thus the crystal analyzer 
placed at ~ 3.4 A"^ {0 = 50°) primarily selects scattering coming 
from region in the sample where p = 2.5 g/cm^. On the other hand, 
different densities are equally weighted in the scattering signal, when 
the crystal analyzer selects ~ 7.4 A"^ at a scattering angle of ^ = 130°. 

Discussion 

The highly charged ions and the resulting strong Coulomb forces 
between them should have further implications for the dynamics of 
the system. As shown in Fig. 3b, the structure factor exhibits a strong 
correlation peak with a maximum of S//(/co) ~ 3.5 at t > 30 ps. 
According to the Hansen-Verlet criterion^^, a first correlation peak 
with Siiiko) > 2.85 indicates a phase transition into a solid with a 
well-defined lattice structure. The measured values in our experi- 
ment clearly fulfill this criterion even if one considers possible errors 
in determining the exact peak value (see Figure 3). Thus, we find the 
situation that we have measured a fluid structure that itself indicates 



that a phase transition into a solid, that is a Coulomb crystal, should 
occur at these conditions. 

Recent models predict that crystalline structures consisting of a 
large number of ions should already have formed in the time span 
between the optical laser and FEL probe^^. Numerical simulations 
give nucleation times of t ~ rca{kBTIM)~^'^ ~ a fewpicoseconds^°'^\ 
where M the ion mass and ~ 10-50 the critical radius for nuc- 
leation (in units of a). Since t is shorter than the characteristic 
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Figure 4 | Ion-ion structure factor calculations. A comparison between 
DFT-MD and MSA predictions for the ion-ion structure factor. In density 
functional simulations both electrons and ions are handled as elementary 
particles. Properties of the electrons are calculated via density functional 
theory using a Mermin functional that accounts for temperature effects 
within the electron subsystem in a statistical sense. The ions instead are 
treated by classical molecular dynamics simulations. This is possible 
because the dynamics of ions and electrons is effectively decoupled with the 
Born-Oppenheimer approximation. By taking snapshots of the ions 
positions, the ion-ion structure factor can be thus calculated. Our ab initio 
calculations are performed with p = 3 g/cm^ and T = 30, 000 K. MSA 
calculations^^ instead use the ionization state (Z) as an additional input. 
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hydrodynamic times (typically tens of picoseconds), the dense 
plasma should have had enough time to evolve into a crystal, a 
process known as Wigner crystallization^^. However, our measure- 
ments show no indication of a crystalline structure which points to a 
state where the ions are locked to a fluid-like structure by strong 
forces. Given the pump-probe difference used, this state exists much 
longer than typical nucleation models predict. 

In order to test the possibility of a crystalline state, we infer 
r > 600 at peak compression with p > 4.5 %/cvs\ and T ~ 
5, 000 K from the hydrodynamics simulations (see Figure 2). At such 
large values of the coupling parameter, a system of ions embedded 
into a screening background of electrons forms a crystalline state 
with a body centered cubic (bcc) lattice^^, as the phase transition is 
found for ~ 400^^. Thus, the highly compressed states in the 
sample are clearly at conditions where we expect a solid. The ratio 
of the scattering power between the solid and liquid phases in the 
{hkl) direction is ^^Psoiid/^^^iiquid ~ i1r(Nd/2)mhki\Fhki\^. Here iVj>l is 
the number of ions in each solid crystallite that has nucleated, f]r is 
the fraction of the sample where T > 400, mh^i and F^ki are the 
multiplicity and the structure factor, respectively, for Bragg reflec- 
tions along the plane of Miller indices (hkl). The hydrodynamics 
simulations (see Figure 2) indicate rjr ~ 8%-10% (for sample densi- 
ties p = 4.5 ± 0.25 g/cm^). Therefore, scattering from this newly 
formed bcc structure should be visible in the data. 

In Figure 3, we have indicated the positions of the Bragg scattering 
peaks for the bcc lattice. For the (112) peak at 4.5 g/cm^ we expect a 
large ratio of the scatter from the solid to the fluid. However, the 
observed scattering at the (112) position remains comparable to the 
signal measured for the (Oil) plane at 2.5 g/cml This indicates that, 
in contradiction to the nucleation models, any significant crystalliza- 
tion has not occurred at the probe time. 

It has been suggested theoretically that the transition to a solid is 
suppressed after rapid changes of the interaction potential between 
the particles leading to a metastable amorphous or glassy state in 
Coulomb systems^ \ One reason for such a behaviour is that the 
transition to a strongly coupled fluid with coordination numbers 
of 3 to 4 is much easier to achieve than the one to a bcc lattice with 
a coordination number of 8. Since we have measured a fluid-like 
structure for conditions for which a bcc lattice should exist, we con- 
clude that we have found compelling indications for such a Coulomb 
glass. Both the Hansen-Verlet criterion^^, which is independent of 
the specific form of the interactions but based on the measured 
height of the first correlation peak, as well as our analysis using the 
radiation-hydrodynamics simulations support this finding. 
Although we can rule out the known solid phases to exist for the 
conditions probed, the determination of the full structure in the 
proposed glassy state will require additional experiments revealing 
the full /c-dependence of the structure factor, as in conventional 
angle-resolved x-ray diffraction. 

The plasma state explored in this experiment (p > 2.5 g/ cm^, and 
T < 1 eV) lies within the parameter range of the envelope of white 
dwarfs with lower luminosity^^ Due to the enhanced ionization state, 
it is also representative of the strongly coupled plasma found in their 
deeper interior, as well as in neutron star matter. The existence of a 
glassy state in white dwarfs would affect their thermal conductivity. 
While the duration of the metastable glassy state is most likely not 
long enough compared to the characteristic cooling time of white 
dwarf stars (which is several billion years), it may nevertheless be 
comparable to the cooling times of neutron stars^^, thus even affect- 
ing their overall luminosity. 
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